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The evolution operator of a discrete-time quantum walk involves a conditional shift in position 
space which entangles the "coin" and position degrees of freedom of the walker. After several steps, 
the coin-position entanglement (CPE) converges to a well defined value which depends on the initial 
state. In this work we provide an analytical method which allows for the exact calculation of the 
asymptotic reduced density operator and the corresponding CPE for a discrete-time quantum walk 
on a two-dimensional lattice. We use the von Neumann entropy of the reduced density operator as 
an entanglement measure. The method is applied to the case of a Hadamard walk for which the 
dependence of the resulting CPE on initial conditions is obtained. Initial states leading to maximum 
or minimum CPE are identified and the relation between the coin or position entanglement present 
in the initial state of the walker and the final level of CPE is discussed. The CPE obtained from 
separable initial states satisfies an additivity property in terms of CPE of the corresponding one- 
dimensional cases. Non-local initial conditions are also considered and we find that the extreme 
case of an initial uniform position distribution leads to the largest CPE variation. 



PACS numbers: 03.67.-a, 03.67.Mn, 03.65.Ud 



I. INTRODUCTION 



The quantum walk (QW) is a reversible process usually 
introduced as a quantum analog of a Markovian process 

Several QW-based algorithms for meaningful prob- 
lems have been developed @, H, H, H, @, 0] and they per- 
form better than the best classical alternatives. In some 
cases, exponential speedups may be obtained [§, [oL [Tot ■ 
There are two versions of QW, based on discrete time 

or continuous time Both have similar dynami- 
cal properties and the later can be obtained by a suitable 
limiting process from the discrete-time walk [IH, [l4| . A 
second level of classification is based on the network over 
which the walk takes place. The QW on a linear chain 
is the simplest possible configuratio n [ill . [TBI . [iHl [ItJ , but 
other topologies such as cycles [H, Il9l]. two-dimensional 
lattices [2(| HH, [23|, [24[ , or n-dimensional hypercubes 

|26| have also been investigated. 



Certain aspects of the QW can be simulated using op- 
tical analogies [13, HH , but to capture the true quantum 
nature of a QW one must deal with entanglement, a key 
resource for quantum information processing. The dis- 
crete time versions of a QW require an auxiliary "coin" 
subspace Tic- The evolution in the position subspace, 
Tip, consists of a conditional shift determined by the 
amplitudes of the states in Tic- The Hilbert space of 
the system is H = He Many QW applications 

involve partial measurements or noise events which se- 
lectively affect a part of the system (frequently, either 
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the coin or position subspaces). The way the system 
is affected depends on the degree of entanglement be- 
tween coin and position just before the event. A precise 
knowledge of this entanglement is required to gain con- 
trol over the long-time behavior of the system, including 
the way it responds to partial measurements or selective 
noise events. It has also been recently suggested that the 
QW protocol may be useful as entanglement generator 
in two-body p9j or even in many-body systems [3C| . 

In a quantum walk the conditional shift operation gen- 
erates entanglement between the coin and position (CPE) 
degrees of freedom of the walker. After several steps, it 
converges to a well defined value which, for a given evolu- 
tion operator, is determined by the initial state [2l|. For 
a one-dimensional (Hadamard) QW, the dependence of 
CPE on the initial state has been characterized using an- 
alytical methods (3l| . However, most algorithmic appli- 
cations require higher dimensions. Furthermore, higher 
dimensional QW's open new possibilities, such as the 
preparation of initially entangled states within the coin 
(CCE) or position subspaces (PPE). 

This paper deals with the characterization of long time 
CPE in a two-dimensional discrete-time QWjSuch as the 
one used in recent algorithmic proposals 0, 0]. This 
system describes the motion of a quantum walker on a 
two-dimensional (2D) lattice or, alternatively, of two in- 
dependent walkers along linear lattices. Our approach 
provides exact results for the two-dimensional QW with 
arbitrary initial conditions and coin operations. We pro- 
vide several examples for a two-dimensional Hadamard 
quantum walk. 

Any real implementation of a quantum system must 
deal with the issue of decoherence, which tends to de- 
stroy quantum correlations. The entanglement decay 
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due to noise has been numerically investigated in one- 
dimensional systems [32| and there are many papers on 
the important subject of decoherent quantum walks, see 
Ref. [33| for a recent review. In this work, we restrict our 
attention to a coherent quantum walk on two dimensions 
and obtain the dependence of the asymptotic CPE on 
the initial coin state for local and non-local initial posi- 
tions. In particular, we show that the non-local case can 
be easily obtained from the local case by adding a weight 
factor to the final integration. 

This work is organized as follows. Section II defines 
the two-dimensional quantum walk and provides the re- 
quired formalism in Fourier space leading to the long- 
time reduced density operator for arbitrary coin oper- 
ations and initial states. In Section III, the method is 
applied to obtain the asymptotic CPE entanglement of 
a two-dimensional Hadamard walk as function of several 
families of initial states. This section includes a discus- 
sion on the additivity of CPE for separable cases. Finally, 
in Section IV, we summarize our results and present our 
conclusions. 



II. ASYMPTOTIC ENTANGLEMENT 

The two-dimensional quantum walk is defined in terms 
of a discrete lattice whose sites are labeled by pairs of in- 
tegers (x,y). One can think in terms of two particles 
moving along two lines or in terms of a single parti- 
cle moving on a plane. For definiteness, in this work 
we adopt the language of a single particle moving on a 
two-dimensional lattice. The set of orthonormal states 
{|a;,y)} spans the position subspace, Tip, of the walker. 
The "coin" degree of freedom is represented by a two- 
qubit space, Tic, spanned by four orthonormal states 
which we label as {\L,L),\L,R),\R,L),\R,R)}. This 
nomenclature is motivated by the quantum walk on a 
line, where \L) and \R) are associated with left or right 
displacements respectively 

The Hilbert space for the system is Ti = Tip (S> Tic ■ A 
generic state is 

4 

i*> = EE£( a 'W)i a? >i'>®tf>. (!) 

®,y 3=1 

where the first sum runs over all lattice sites and we used 
a compact notation j = 1 . . .4) for the coin states 
{\L,L),\L,R),\R,L),\R,R}}. One step of the evolution 
is described by 

|¥(n + l)> = tf|*(n)>, (2) 

where n is a step counter and the evolution operator is 

U = S-(I P <E> U c ) (3) 

with Ip the identity operator in Tip. The evolution com- 
bines a unitary coin operation Uc in Tic with a shift 



operator 

S = - l >v)fav\ ® + \x,y + l){x,y\ ® |2)(2| 

x,y 

+\x,y-l)(x,y\ ® |3)(3| + \x + l,y)(x,y\ ® |4)(4|} (4) 

which performs the conditional displacements deter- 
mined by the coin state. The correspondence between 
coin states and displacements is not unique. The shift 
operator defined in eq. (0| represents the extension of a 
one-dimensional walk to two dimensions in a 45° rotated 
lattice. This choice for S simplifies the description in the 
Fourier representation, without loss of generality. 

Each spatial component of the wavevector, \ip x ,y) = 
(x,y\V) = J2 4 j= i fj(x,y)\j), evolves as 

|^»,»(n + l)) = Mi|^»+i lW (n)) + M 2 |^x, v - 1 (n)) (5) 
+ M 3 \ip XtV+ i(n)} + M 4 \ip x _ ltV (n)}. 

in terms of four operators (Mi) acting in TLc- For a 
generic coin operation, Uc — Ci.jK)OI; these oper- 
ators are 

4 

M i = \i)(i\U c = J2 C ^\ i )(j\- 

7=1 

The Fourier transform, as first noted in this context 
by Nayak and Vishwanath [ll[ , is extremely useful when 
single-step displacements are involved because the evolu- 
tion operator is diagonal in fc-space. The Fourier trans- 
form of the position eigenstates is 

|k> = (7) 

r 

where r is the vector with integer components (x, y) and 
k is a vector with real components (k x , k y ) in the interval 
[— 7r, 7r]. The k th component of the wavevector, eq. (TjQ), is 

4 

m^fjo^m (8) 

3=1 

has amplitudes fj(k) = Yl r e ~ r fj( r )- ^ s mentioned 
before, the linear map in eq. ([5]) can be expressed as the 
action of a diagonal operator in k-space 

|Vfc(n + l)) = U k \^ k (n)) 

= (e-' fc "M 1 + e^M 2 (9) 
+ e- ik yM 3 + e ik *Mt) |Vfc(»)>. 

This equation defines a unitary operator U k which is rep- 
resented by a 4 x 4 matrix. The basic idea behind our 
approach is to use the spectral decomposition of U k to 
obtain information about the long-time evolution of the 
system. 
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Let us consider the eigcnproblcm for the unitary oper- 
ator U k with eigenvalues e lu)k and corresponding normal- 
ized eigenvectors 

4 

|eH))=£ aj -(w fc )|j). (10) 

3=1 

Using the spectral decomposition for the state of the 
system after n steps is 

\Mn)) = EWfc(O)) = £ ^"^H) l e M>> (11) 
where the sum is over the set of eigenvalues of U k and 

4 

F( Wfc ) = (e(o; fe )|^(0)) =E"K^)/ J -(k). (12) 

i=i 

Note that /j(k) are the Fourier-transformed initial am- 
plitudes, so F(uj k ) contains all the information about the 
initial state. 

Since we deal with pure states only, we use the von 
Neumann entropy of the reduced density operator, or 
entropy of entanglement, defined as 

E = -trace (p c log 2 p c ) (13) 

as a measure of coin-position entanglement (CPE). In 
this expression, p c — tracep (p) is the reduced density 
operator obtained from p = U n poU jt n by tracing out the 
position degrees of freedom. Since p c has dimension 4, 
this quantity is E £ [0, 2], i.e. E = for a product state 
and E = 2, for a maximally entangled state. For a quan- 
tum walk on a line, p c is two-dimensional and E € [0, 1] 
and for p c with dimension d, E is in the interval [0, log 2 d] . 

The inverse Fourier transform required to generate p at 
arbitrary times can not be computed exactly, but since 
our main interest is to obtain the entropy of entangle- 
ment, we can avoid this issue by using Parseval's theo- 
rem. The reduced density operator, p c , is obtained by 
taking the trace in fc-space of the density operator 

p c = tracer (p) = I — \i>k) {$k\, (14) 

where tracer (•) traces over (k x ,k y ) and the integration 
in rf 2 k = dk y dk x has limits [— n, n}. This expression can 
be evaluated after many steps of the evolution in a form 
completely analogous to the one used in Ref. [3l[ for a 
quantum walk on the line. From cq. pi[) . after n steps, 

KW fe |= £ e*"-~ u '^ n F{u k )F*{J k )\e{u h )){e{u' k )\. 
Wk-^' k } 

In the asymptotic limit n> 1, according to the station- 
ary phase theorem, only terms with tu k = u>' k contribute 
in eq. (fl"4|) as discussed in detail in Ref. [ll| . Thus 

k = /S £ \F(u k )\ 2 \e(u; k )){e(u k )\, (15) 



where we use a caret (") to indicate that the asymptotic 
limit p c = limn^oo p c (n) has been taken. According to 
eq. ()10|) . the matrix elements of p c are 

(i\Pc\j)= J S^(k) (16) 

in terms of 

VijQa) = £ I^K)| 2 a 4 (w fc )a*K). (17) 

Note that these expressions satisfy Vi,j = ~P*j, as re- 
quired by the hcrmiticity of p c . 

At this point we need to specify the form of the initial 
state. Since the main interest of this work is to char- 
acterize the long-time coin-position entanglemet (CPE) 
generated by the evolution of the quantum walk, we shall 
consider only separable coin-position initial states, i.e. 

|*(0)> = W) ® |x) (18) 

with initial position \ip) = X) r a ( r )l r ) an< ^ initial coin 
|x) = J2j c j\j)- These states have no CPE and 
fj{r) = Cj a(r) in cq. (fTJ) , so that their ^-component is 

|V fc (o)) = a(k) | X ) (is) 

in terms of the Fourier-transformed initial amplitudes 
a(k) = Y^r e ~ Ta ( r )- F rom ec l- l|12p one readily obtains 
for this kind of initial conditions, 

|FK)| 2 = |&(k)| 2 £a> fe )a ; K) Cj cr. (20) 

3,1 

From this expression for a given coin operation and ini- 
tial state one can obtain p c and, after diagonalization, 
the entropy of entanglement E, from eq. (|13[) . Thus, the 
dependence on the initial conditions of the CPE gener- 
ated by the evolution after many steps can be explored 
using this method. We emphasize that it can be gener- 
alized in a straightforward form to arbitrary dimensions 
and coin operations Uc, as long as one is able to solve 
the eigenproblem for U k analytically. 

III. DEPENDENCE ON THE INITIAL STATE 

In order to discuss the dependence of CPE on 
the initial state, we must specify a coin opera- 
tion. We shall consider the case of a two-dimensional 
Hadamard walk. When applied to the basis states 
of a single qubit, a Hadamard operation generates 
the balanced superpositions H\L) = (\L) + \R))/y/2 and 
H\R) = (\L) — \R))/^/2. This coin operation is a com- 
mon choice in the literature of one-dimensional quantum 
walks [ill, EH 113 and information on the asymptotic 
coin-position entanglement (CPE) level is available for 
this case [U [3l| . Our method is also applicable to other 
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cases of interest, such as the Grover or DFT coin opera- 
tions [2l|. 

For a quantum walk in two spatial dimensions the coin 
operation Uq = H ® H is a natural extension of the 
Hadamard walk on a line. In this Section, we apply the 
general formalism previously outlined in the previous sec- 
tion to this case and characterize the asymptotic CPE 
level for several initial conditions. 

From eqs. (|6]) and ([9]), the explicit form for the operator 
U k forU c = H®H is 



U, = - 







e -iK 


g ik x 








gifcy 


~—ik y 


g— iky 


— p~^v 


— e~' lk 


e ik x 


pik x 


pik x 





\ 

) 



(21) 



The eigenproblcm for this operator can be solved exactly. 

, e _lw - } with 



Its four eigenvalues are {e luJ+ , e e 1 ^- 



COSOJ-I- 



cos k x — cos k y ± v Afe 



(22) 



Afe = cos 2 k x + 6 cos k x cos k y + cos 2 k y + 8. 

As usual, we shall refer to the phase of any of these eigen- 
values by the symbol Wfc. The normalized eigenvectors 
are of the form (p~0|) with components 



1 

N 
1 

N 
1 

N 



:\l-e 2luJk ] (23) 

— 1 + e iuJk (e^ 1 — e ifCy ^ + e 2iuJk e^^^v) 

— 1 + e iuJk (e^ 1 — e~~ ifcy ) + e 2iuJk e^ fc:r— ^v) 



no 
0.3 
04 



— [l + e 1 ^ (e lfc » + e _tt ») + 



The normalization constant N is a positive real, chosen 
so that Y^=i \ a i\ 2 = 1- 



A. Separable initial states 

Let us first consider in the detail the simple case of 
an initial state localized at the origin (x = y = 0) with 
initial coin state \x) = \L,L), or c\ = 1 and Cj = for 
j 7^ 1. For a localized state, a(k) = 1 and the projection 
on k-space, eq. (fl"8j) . is simply = Thus in this 
case, eq. (|20|) reduces to 



|F(w fc )| 2 = 2(1 - cos(2w fc )). (24) 
This expression is used in eq. (|17p to find 

(25) 

A^" 1 (cos 2 k x + 4 cos fej, cos fc^ + cos 2 fcj, + 3) . 



P M (k) = V|F(c fe )| 2 -|«iK)| 2 



After averaging over k, the matrix clement of the reduced 
density operator is obtained, 



p c (i,i) = J 0P w (k)EC 1 = 



(26) 



The other independent elements of p c may be calculated 
in the same form, 

5 - 3V2 



Mi, 2) 


= /Sc(l,3) = 


c 2 


Pc(l,4) 


= /5 c (2,3) = 


c 3 


/5 C (2,2) 


= /5 c (3,3) = 


c 4 


Pc(2,4) 


= Pc(3,4) = 


c 5 



3-2^ 
8 

2V2- 1 
8 

V2 - 1 



(27) 



The eigenvalues for p c are Xi = 1/2, A 2 = 4C3 and 
A 3 = A 4 = 4C5, so that the asymptotic entropy of en- 
tanglement for this case is 

E = -trace (p c log 2 p c ) ~ 1.744 

This quantity can be compared to the asymptotic entan- 
glement of a one dimensional Hadamard walk that starts 
at the origin with \L) (or \R)) as the initial coin. In this 
cases, the asymptotic CPE is E ~ 0.872, [H|,[3l||. Tnus , 
for \L,L) wc obtain exactly twice CPE as in the one- 
dimensional case. This is due to the fact that both the 
initial coin state and the coin operation H ® H are sep- 
arable and the quantum walk on the plane decomposes 
into two independent one-dimensional walks. In these 
cases, an additivity property applies for the asymptotic 
CPE as we discuss in more detail below. 

Let us suppose that the initial state is separable with 
respect to both walkers, i.e. 



I*) 



(28) 



where = |^) <g> for i = 1,2. The states \cf>i), 
spanned by {|x}}, describe initial positions and the states 
\Xi), spanned by {\L), \R)}, are one-qubit initial coin 
states. If the coin operation can be written as a product 
two one-qubit operators, Uc — Uci <%> Uc2, the evolution 
operator U defined in eq. ([3]), is separable, U = U x ®U y , 
with U x — S x ■ (I x ® Uci) and a similar expression for 
U y . Then, the separability of the initial state is preserved 
by the evolution. After n steps, the state vector can be 
expressed as 

|*) = t^|*i(0))®L^|* a (0)), (29) 

and the reduced density operator 

Pc = traccp (l^) (^|) = pci ® PC2i (30) 

is also separable, with pci = trace x (l^i) |) and a 
similar expression for pci- The subadditivity property 
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of the von Neumann entropy for a separable density op- 
erator [3f| implies, 



E(p C i ® Pci) = E(pci) + E(p C2 )- 



(31) 



In other words, for a separable initial state and separa- 
ble coin operation, the CPE for a two-dimensional QW 
can be obtained from the CPE of the underlying one- 
dimensional quantum walks. 

One can use this property as a witness of correctness 
of our method for the two-dimensional case. Let us con- 
sider a a localized initial position with a generic separable 
initial coin state 

|*(0)) = |0, 0) ® | X i(0i, <h)) ® |xa(02, fc)) (32) 

where = cos6j \L) + e l( ^ j sindj \R) for j = 1,2 are 
generic one-qubit states. The four real parameters that 
specify this initial state are restricted to the intervals 
0j G [-7r/2,7r/2] and fa G [-7T,7t]. 

According to eq. (|3~Tj) . the asymptotic CPE in 
this case can be computed from the CPE of a 
one dimensional quantum walk with initial condition 
|*i(0)) = |0) (g> \xi{0,4>)). As mentioned before, for this 
case an exact expression for the CPE has been obtained 
in Ref. [3l| . in terms of the initial coin state and the 
quantity 



A(M) 



1 + (1 - 4(A - 2b\ sin 49 cos </>)) ; 



(33) 



The constants in this expression are Ao = (y/2 — l)/2 
and bi = (2 — y/2)/4 and the asymptotic CPE for this 
one-dimensional, localized case is given by (3lj 



-Alog 2 A-(l-A)log 2 (l-A). (34) 



Using additivity, one can obtain the CPE for localized, 
separable initial conditions (eq. I32p from this expression 
as 

E(0 1 ,<f> 1 ;0 i ,<h) = E 1 (6 1 ,<f> 1 ) + E 1 (e2,(h). (35) 



The RHS of this equation can be evaluated using eqs. 
and (|33[) and the LHS can be evaluated as described in 
the last section. Thus, eq. (|33|) . which applies to sepa- 
rable localized initial conditions only, is a useful witness 
for the correctness of the proposed method. 

We shall now discuss some specific examples of the 
CPE dependence on the initial coin state. Let us consider 
a particular instance of eq. (|32[) with 9\ = and rename 
(#2,^2) — > (0,0), so that the initial separable coin state 
is 

| X (/) (0,0)) = \L) ® (cos6\L) +e^sm6\R)) . (36) 

For these states a simple calculation based on eq. (f2"0)) 
leads to 

\F\ 2 = |ai| 2 cos 2 6»+ |a 2 | 2 sin 2 6» + sin(26»)K(ata 2 e-^) . 
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Figure 1: (a) Asymptotic entropy of entanglement (CPE), 
E{9,(f>), for a localized separable initial state, eq. (|32|l . with 
the particular initial coin state defined in eq. (|36[) . (b) Con- 
tour plot of the same surface. 



This expression is used in eq. (|T7|) to evaluate Vij (k) . Af- 
ter averaging in k-space, the dependence of the asymp- 
totic reduced density operator on the initial state param- 
eters can be expressed in terms of the functions 

h(6,4>) = sin(26») cos 

f{9,4>) = h{9, 4>)+ cos(20) + 1 (37) 
g(6,(p) = f(9,(/>)-iV2sm<t>sm(26)-l, 

and the constants Ci defined in eqs. (f2"o| - (f2"T)) . The ex- 
plicit form of the reduced density matrix is 

C\ + C 2 f C 2 g C 5 + C 3 f C 3 g 

C 2 g* Ci-C 2 f C 3 g* C 2 -C 3 f 
1 C 5 +C 3 f C 3 g i + C 5 / C 5 g 

C 3 g* C 2 -C 3 f C 5 g* C A ~C b f. 
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-71/2 -71/4 TC/4 TC/2 



1 

Figure 2: Contour plot for the CPE, i?(#i,#2), obtained from 
the initial separable state eq. (|32[) , with <j>i = <f)2 ~ 0. 

The exact diagonalization of this operator leads to the 
asymptotic entanglement as a function of the initial pa- 
rameters, E(9,<fi). Fig. [T] shows the entropy of entangle- 
ment as a function of initial coin state defined in eq. (|36p . 
Alternatively, the surface shown in Fig.[T]can also be cal- 
culated from the RHS of eq. (|35|) . 

E(e,4>)=E + E 1 (6,ct>) (38) 

with Ei given by cq. (34|) and E Q = E^O, <j>) ~ 0.872. 

In Fig. [2] we show the CPE for the case in which both 
initial coin states are superposition states. We have set 
cf>i = fa = in cq. (|32|) and calculated E(6i, 9 2) following 
the method outlined in the previous section. The result 
is consistent with the additivity property, eq. (|35p . 

The previous argument, based on the separability of 
the motion, may be extended to any number N > 2 of 
quantum walkers with a separable initial state, |^(0)) = 
|*i(0)) ® . . . ® |*jv(0)), evolving under a separable coin 
operation, Uc = A\®. . .®An. Then, the CPE can be ob- 
tained by addition of the corresponding one-dimensional 

CPE's, £"(*(())) = J2f=i E i(^i(°))- The maximum en- 
tanglement of E = N can be obtained of all initial states 
are properly prepared. 

B. Entangled initial states 

Up to this point we have discussed separable initial 
states only. Let us now discuss the effect on CPE of ini- 
tially entangled states. Strict subadditivity of the von 
Neumann entropy holds only for separable density op- 
erators, so in this case CPE can not be obtained from 
one-dimensional calculations. Let us start with localized 
states 

|*(0)) = |00) ® Ix) (39) 



and consider two families of initial coin states, 

|x (//) ) = cos 9 \L, R) +e*+ sin 9 \R,L) (40) 
\X {IH) ) = cos6\L,L)+e i4 'sme\R,R) (41) 

which describe entangled states in Tic- Their entropy of 
entanglement is defined as 

S = -tracer [(|x>0d)log a (|x><xl)] , (42) 

with the partial trace taken over either of the one-qubit 
coin subspaces spanned by {\L), \R)}- This quantity is 
normalized to one and measures the initial coin-coin en- 
tanglement (CCE) in |x). As shown in Fig. 0] (b), it 
depends on the parameter 9 alone. Note that, since the 
coin operation Uc — H ®H is separable, the initial CCE 
is preserved by the evolution. In particular, the maxi- 
mally entangled Bell states 

|*±) = -L(\L ) R)±\R,L)) (43) 

I**) = ±(\L,L)±\R,R)), (44) 

included in these families, appear associated with maxi- 
mum or minimum asymptotic CPE values. 

For the first family, |x^^)j straightforward evaluation 
of eq. (HOI) leads to 

\F\ 2 = \a 2 \ 2 cos 2 9 + \a 3 \ 2 sin 2 9 + sin(20) 3? [a* 2 a 3 e'^] . 

A similar expression holds for Ix^ 111 ^) with (012,0(3) re- 
placed by (a\, 0(4). These expressions are used in eq. (fT7|> 
to obtain the elements of the hermitic matrix "P(k). Af- 
ter the k-average is done, the long-time reduced density 
operator p c is obtained. For Ix 1,11 ^), Pc has seven inde- 
pendent elements which can be expressed in terms of the 
constants C, and the functions h, f, g defined in cqs. (|37p . 
as 

p c (l,l) = Pc(4,4) = d + C 3 h 
p c (2,2) = C 3 (f-2h) + ±(f-h+l) 

p c (3,3) = -Csf + ^ih-fi+d (45) 

p c (l,2) = -p* c (2,4) = C 3 (g* + l)-C 2 (f-h) + C 5 

p c (l,3) = -p* c (3,4) = C 3 (g+l) + C 5 (f-h)-C 2 

p c (l,4) = -C 3 (h + l) 

6,(2,3) = [3(fr-l) + 2(/ -<?*)]. 

Similar expressions can be obtained for the other entan- 
gled family, |x^ //7 ^)- After diagonalization of p c , an ex- 
act expression for the asymptotic CPE E is obtained. 
For both cases, the dependence of this quantity on the 
initial coin state parameters {9, <f>) is shown in Fig. [3l It 
varies in the (approximate) range [1.744, 1.978], with its 
minimum associated to separable initial coin states and 
its maximum associated to fully entangled initial coin 
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e 



Figure 3: Asymptotic entropy of entanglement resulting from 
entangled initial coin states. Upper panel: initial coin in fam- 
ily |x' /J ^) defined in eq. (|40|l . Lower panel: initial coin in 
family |x ( " 7) ), defined in eq. gTJ. 



states. Asymptotic CPE and initial CCE are related in 
the sense that maximum CPE is associated with maxi- 
mum CCE while minimum CPE is associated with initial 
product states (i.e. no CCE). However, maximum initial 
CCE does not imply maximum asymptotic CPE; for the 
family x^ /7 \ maximum CPE is obtained from the Bell 
state I*" 1 "), but l^ - ), also fully entangled, leads to the 
intermediate CPE value E ~ 1.888, see the upper panel 
Fig. [3J For the symmetric family X i the same val- 
ues are obtained but in this case |$~) yields maximum 
entanglement and \& + ) leads to an intermediate value 
E ~ 1.888, see the lower panel Fig. H 

A comparison between both panels in Fig. [3] shows 
that the relation E{x {III) ;0,<l>) = E{x (II) ]-0,(f) is sat- 
isfied. From a mathematical point of view this may be 
traced to the fact that the reduced density operators 
Pc{x [III) \&A) and p c (x (//) ;-6»,0) have the same eigen- 
values. We also note that both initial coins have the same 
CCE, as indicated by S(9), which is an even function (see 
Fig. ffl b). 



1.978 

1.888 
E 
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-71/2 -7C/4 7t/4 TC/2 
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_ (b) 






CCE 











-7t/2 -tc/4 Tt/4 tc/2 



e 

Figure 4: (a) Asymptotic CPE, E{0,<j>), for the families of 
initial coin states |x ) ( re( h dashed line) and |x } (blue, 
full line), (b) Initial CCE S{9), defined in eq. (42j), for the 
same families of states. 

Up to this point, only localized initial conditions have 
been considered. In the one-dimensional case it is known 
that non-local initial states can modify significantly the 
dynamics of the QW. For instance, the Survival Prob- 
ability (the probability of finding the walker in a given 
region which includes the starting point) decays as t~ l 
for localized initial states and as t~ 3 for non-local initial 
states with appropriate relative phases [36[ , so the prob- 
ability flux spreads out faster than in the localized case. 
In the specific case of asymptotic CPE, previous work 
for a one-dimensional QW suggests that non-local initial 
states can give rise to a broader range of variation for 
CPE levels [3l|. However, this initial work was limited 
in scope because the one-dimensional case does not allow 
for initially entangled (coin or position) states as the 2D 
case does. We now consider initial sates of the form 

|*(0)) = |V(«,j3)>® \X) (46) 

where the coin state |x) is chosen within the family 
Ix^') as either one of: the Bell states |x) = | V I >± ) de- 
fined in eq. (|43[) or the product state \LR). The initial 
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Figure 5: Asymptotic entropy of entanglement E for the non- 
local separable initial condition ^ s (a,/3), eq. (|47[1 . with non- 
locality restricted to the £- direct ion. The surfaces correspond 
to different initial coin states: (a) |^ + ), (b) \LR) and (c) 
l^ - ), defined in eq. (|43[) . In each case, the flat surface indi- 
cates the CPE level for a localized initial position with the 
same initial coin state. 

position \ip(a,(3)} is a superposition state on the two- 
dimensional lattice. We start by considering the effects 
of small non-locality, with the initial states 

\ip s ) = (cosa| - l) x + e ?;,3 sina|l) 2; ) ® \0) y (47) 
\ip e ) = cosal - 1,1} +e i/3 sina|l,-l) (48) 

which describe non-local separable (jp s ) and entangled 




Figure 6: Asymptotic entropy of entanglement E for the non- 
local entangled initial condition tp e (a,f3), eq. (|48|l . The sur- 
faces correspond to three different initial coin states as de- 
scribed in Fig. [S] In each case, the flat surface indicates the 
CPE level for a localized initial position with the same initial 
coin state. 



(ip e ) states in Hp. In the later case, the degree of bi- 
partite entanglement (PPE) is given by the entropy of 
entanglement, S(a), which can be read from Fig. [4] (b) 
with the replacement 6 — > a. In the next subsection we 
shall consider the effect of extended non-locality on CPE. 

For these non-local initial conditions the Fourier trans- 
formed amplitudes, a(k) ~ X)r e_ik rfl ( r )' mu st be cal- 
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culated. The relevant expressions are 

|a s (k)| 2 = 1 + sin(2a) cos(2/c K 
|a e (k)| 2 = 1 + sin(2a) cos[2(fc x 



0) (49) 

-ky)+0\. 



The calculation of the reduced density operator is sim- 
plified by noting that eqs. |[T7]) and (|2TJ)) imply 



V'(k) = |5(k)| 2 P(k). 



(50) 



Thus, the matrix T"(k) for the non-local case is expressed 
in terms of the corresponding matrix "P(k) for the local 
case with the same initial coin state. The correspond- 
ing density operators are obtained from integration in 
k-space of the 'P(k) previously calculated for the local 
case, with integrating factors given by eqs. (|49|) . 



Pc 



/g|5(k)|^(k). 



(51) 



Diagonalization of this operator leads to the correspond- 
ing asymptotic entropy of entanglement. E(x',cx,(3), 
which depends on the initial coin state and the two pa- 
rameters that specify the initial position. 

A comparison of the resulting CPE surfaces E(a, /3) 
with the corresponding localized cases values of E show, 
in most cases, variations due to position non-locality. 
The results for separable non- locality ("0 S ) are shown in 
Fig.[5]and those for entangled non-locality, (ip e ), in Fig. [5] 
The flat surfaces represent the CPE of the localized cases 
with the same initial coin states. The largest variations 
due to non-locality are associated to balanced superpo- 
sitions, a = ±7r/4. Notice that the relative phase of the 
initial coin state can cause significant changes in CPE. 
The overall effect of this restricted non-locality on CPE 
is small, but this was to be expected since we are consid- 
ering only small derealizations about the origin. In the 
next section we consider the effect on CPE of non-local 
states with a large spread over the 2D lattice. 



C. Extended non-locality 

We have considered initial position states with a small 
degree of non-locality, i.e. the initial amplitudes that 
are non-zero only in a few sites around the origin and 
found that the changes in asymptotic CPE due to ini- 
tial non-zero amplitudes at sites i = ±1 or 1/ = ±1 
as compared with the local case x = y = 0, are small. 
However, is more extended non-locality capable of sig- 
nificant changes on the long-time CPE? Let us approach 
this question, considering an isotropic Gaussian position 
distribution for the amplitudes, a(r) oc e~( x +y ^ 2a , 
where a > is a characteristic width. The Fourier- 
transformed amplitudes arc 5(k) = Ce" 5 ^ 1 ^'" with 
C a constant obtained from the normalization condition 
J^|a(k)| 2 = l. Since 

lim e~i^+ k > 2 = ^-6(k x )6(k y ) 
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Figure 7: Upper panel: contour plot for E(9 
X (0, 4>) defined in eq. (|40|l and initial position uniformly 
distributed over the plane. Lower panel: constant cj> sections 
for the same surface. 



where S(-) is a Dirac's delta function, extended non- 
locality in position transforms the integration factor 
|a(k)| 2 in a product of two delta functions and, in this 
limit, the integration in eq. (|5ip becomes trivial. The 
resulting reduced density operator is 



(52) 



Thus for the extended position distribution, the 
asymptotic CPE depends only on the initial coin state 
which determines |-F| 2 . Let us choose the initial coin 
state within the family of entangled states x^K^'Wj 
defined by eq. (|4"0l . In this case, the elements of the 
long-time reduced density operator, p c , result from eval- 
uating the expressions for V{k) obtained in the localized 
case, for k x = k y = 0. The result can be expressed in 
simple form. The diagonal elements arc 



p c (l,l) = (sin(20)cos<73 + 3)/16 

p c (2,2) = (1 - sin(26>) cos </> + 8 cos 2 6l)/T6 

p c (3,3) = (9-sin(26»)cos</>-8cos 2 6»)/16 

/5 C (4,4) = (3 + sin(26»)cos^)/16, 
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with trace (p c ) = 1, as expected. The expressions for the 
off-diagonal elements are, 

pcO-,2) = (l-4cos 2 6» + sin(26i)cos0)/16 

/5 C (1,3) = (-3 + sin(2<9)cos0 + 4cos 2 6>)/16 

Pc(lA) = (-1 + sin(26»)cos0)/16 

Pc(2,3) = (-1 + sin(26»)cos0)/16 

ft. (2, 4) = (-1 -sin(26»)cos0 + 4cos 2 6»)/16 

p c (3,4) = (3-sin(26»)cos</>-4cos 2 6')/16. 

The eigenvalues of this operator can be obtained analyt- 
ically and the entropy of entanglement calculated as a 
function of the initial coin state. The resulting surface, 
E(6, </>), is shown in the upper panel of Fig. [7J Two as- 
pects of this surface are remarkable: (i) full asymptotic 
entanglement E max = 2 results for the initial coin |^ + ). 
Note that this maximum is basically flat in the <fr direc- 
tion, so it is robust against small variations in the relative 
phase of the initial coin state. And (ii) the minimum en- 
tanglement is now (exactly) E m i„ — 1 for the initial coin 
l^ - ), as shown in the lower panel of Fig. [7J Finally, the 
initial product state \L,R) results in a low CPE level 
of E ~ 1.20. Comparison with the local case in Fig. [3] 
shows that extended non-locality increases significantly 
the range of variation for asymptotic CPE. This result 
also contextualizes the small variations in CPE due to 
initial non-local positions in the neighborhood of the ori- 
gin, as the variations may be larger when more extended 
initial states in position space are considered. 

IV. SUMMARY AND CONCLUSIONS 

This work describes a method for the exact character- 
ization of the long-time (asymptotic) coin-position en- 
tanglement (CPE) of a discrete-time quantum walk on 
a two-dimensional (2D) lattice or, alternatively, of two 
independent walkers on a line. In order to quantify the 
bipartite entanglement of the pure state p, the von Neu- 
mann entropy of the reduced density operator, p c is used. 
This quantity, E(p c ), is scaled so that it varies between 
for a product coin-position state to 2 for a fully en- 
tangled state. This corresponds to the [0, 1] variation 
of the one-dimensional case. The initial condition is all- 
ways chosen as a coin-position product state, so that no 
CPE is initially present. With this sole restriction, the 
general formalism leading to the exact calculation of the 
long time reduced density operator is presented. This 
treatment allows both local and non-local initial posi- 
tions and arbitrary coin operations. The expressions we 
have presented can be readily applied to quantify the 
CPE of quantum walks in higher dimensions, such as the 
n-dimensional hypercube pa . l2q |. The exact nature of 
our results allow us to distinguish between small vari- 
ations in CPE and, for instance, clearly identify initial 
states which lead to full entanglement. 

In order to illustrate the kind of results that can be 
obtained, we have considered in detail the case of a 



Hadamard coin operation, Uc — H ® H . Similar cal- 
culations can be done for any coin operation (i.e. Grover 
or DFT coins) for which the relevant eigenproblem in 
k-space has been solved. We have first considered the 
case of localized positions with separable initial coin 
states. Then each walker starts in well defined states 
in the coin and position subspaces and the motion under 
a Hadamard coin remains separable. In this case, the 
CPE satisfies an additivity property which allows it to 
be expressed as the sum of the CPE for the correspond- 
ing one-dimensional motions. This additivity property 
can be generalized to the case of N independent walkers 
with separable initial condition and coin operation. The 
CPE obtained is rather high, above 73% of the maximum 
value in all cases considered, except for the extended non- 
locality discussed below. The maximum CPE depends 
on wether one can tune both initial coins or not. If both 
coins are tunable, the maximum CPE, E = 2, can be 
obtained. 

When entangled initial coin or position states are con- 
sidered, the motion is no longer separable and the calcu- 
lation of CPE can not be reduced to the one-dimensional 
case. In this case, there is the (coin-coin or CCE) en- 
tanglement present in the initial state and the long-time 
CPE generated by the evolution. These are new results, 
as previous analytical work based on one-dimensional 
walks [3l| did not allow for this possibility. Extreme 
values for CPE (maxima or minima) appear associated 
with maximum values of CCE. Asymptotic CPE and ini- 
tial CCE are related in the sense that maximum CPE 
implies a maximally entangled (CCE) initial coin state 
and minimum CPE implies an initial product state (i.e. 
no CCE). However, there arc maximally entangled initial 
coins which do not lead to maximum CPE. 

We have also consider the effect of non-local initial 
conditions. We have shown that the calculation of the 
asymptotic density operator for this case can be reduced 
to that of the local case with the same initial coin state, 
provided an additional integration factor is included in 
the final step. Initial superpositions of two close sites 
can produce small changes in CPE (smaller than 1%) 
with respect to the localized case. However, when ex- 
tended non-local states are considered, large variations of 
CPE with respect to the local case may be obtained. We 
have investigated an initial state with Gaussian ampli- 
tude distribution about the origin. In the uniform limit, 
the variation for CPE is large. Values between E = 1 
and E = 2 can be obtained depending on the initial coin 
state. 

Many studies of entanglement are restricted to the few- 
qubit case. The 2D quantum walk involves a two-qubit 
coin subspace and a position register which must have 
several qubits in order to accommodate the evolution 
during a significant number of steps. Thus, the entan- 
glement in these systems is a complex issue, involving 
different kinds of bipartite entanglement (some of which 
we code-named as CCE, PPE or CPE). The analytical 
results presented here are a step towards a better under- 



11 



standing of the dynamics of entanglement in quantum 
walks. The relation between asymptotic CPE and the 
initial state is important, for instance, in order to gain 
insight on how a partial measurement will affect the walk. 
It also allows one to choose the initial conditions that will 
lead to the desired entanglement level. Entanglement is a 
key resource for quantum information processing, so we 
hope that these and similar results may stimulate new 
algorithmic applications. 
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